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ABSTRACT 

Context. PAHs are thought to be a ubiquitous and important dust component of the interstellar medium. However, the effects of their 
immersion in a hot (post-shock) gas have never before been fully investigated. 

Aims. We study the effects of energetic ion and electron collisions on PAHs in the hot post-shock gas behind interstellar shock waves. 
Methods. We calculate the ion-PAH and electron-PAH nuclear and electronic interactions, above the carbon atom loss threshold, in 
Hn regions and in the hot post-shock gas for temperatures ranging from 10 3 - 10 8 K. 

Results. PAH destruction is dominated by He collisions at low temperatures (T < 3xl0 4 K), and by electron collisions at higher 
temperatures. Smaller PAHs are destroyed faster for T < 10 6 K, but the destruction rates are roughly the same for all PAHs at higher 
temperatures. The PAH lifetime in a tenuous hot gas (n H a 0.01 cm -3 , T a 10 7 K), typical of the coronal gas in galactic outflows, is 
found to be about thousand years, orders of magnitude shorter than the typical lifetime of such objects. 

Conclusions. In a hot gas, PAHs are principally destroyed by electron collisions and not by the absorption of X-ray photons from the 
hot gas. The resulting erosion of PAHs occurs via C2 loss from the periphery of the molecule, thus preserving the aromatic structure. 
The observation of PAH emission from a million degree, or more, gas is only possible if the emitting PAHs are ablated from dense, 
entrained clumps that have not yet been exposed to the full effect of the hot gas. 
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1. Introduction 

The mid-infrared spectral energy distribution of the general in- 
terstellar medium of galaxies is dominated by strong and broad 
emission features at 3.3, 6.2, 7.7 and 11.3 //m. These fea- 
tures are now univocally attributed to vibrational fluorescence 
of UV pumped, large 50 C-atoms) Polycyclic Aromatic 
Hydrocarbon (PAHs) molecules. These large molecules are very 
abundant (3 x 10~ 7 by number relative t o H-nuclei) an d ubiqui- 
tous in the ISM (for a recent review see [Helens 2008). Besides 
large PAH molecules, the spectra also reveal evidence for clus- 
ters of PAHs - containing some hundreds of carbon atoms - and 
very small grains (^ 30 A). Indeed, PAHs seem to represent the 
extension of the i nterstellar dust size distribution into th e molec- 
ular domain (e.g. lDesert et al.|["l990HDraine & H I2001 . 

PAH molecules are an important component of the ISM, for 
example, dominating the photoelectric heating of neutral atomic 
gas and the ionization balance of molecular clouds. Small dust 
grains and PAHs can also be important a gents in coo ling a hot 
gas, at temperatures above ~ 10 6 K (e.g. Dwekl[T98~7l) . through 
their interactions with thermal electrons and ions. The energy 
transferred in electron and ion collisions with the dust is radiated 
as infrared photons. The evolution of dust in such hot gas (T £ 
10 6 K), e.g., within supernova remnants and galactic outflows, is 
critical in determining the dust emission from these regions and 
therefore the cooling of the hot gas. The destruction of PAHs and 
small dust grains in a hot gas may also be an important process 



in the lifecycle of such species (IDwek et al.lll996t iJones et al.l 
fl996h . 

Observationally, there is little direct evidence for PAH 
emission uneq uivocally connected to the hot gas in super- 
nova remnants. iReach et all (120061) have identified four super- 
nova rermTants_with_JR colors that may indicate PAH emis- 
sion. iTappe et ail ((2006) have detected spectral structure in the 
emission characteristics of the supernova remnant N132D in the 
Large Magellanic Cloud that they attribute to spectral features 
of PAHs with sizes of 4000 C-atoms. Bright 8 fim emission 
has been observed by IRAC/Spitzer associated with the X-ray 
emission from the stellar winds o f the ionizing stars in the M 
17 Hn region dPovich et al.ll2007T) . Likely, this emission is due 
to PAHs - probably, in entrained gas ablated from the molecu- 
lar clouds to the North and West of the stellar cluster. Finally, 
bright PAH emission has been detected associated with the hot 
gas of the galactic wind driven by the starburst in the nucleus 
of the nearby irregular galaxy, M 82 dEngelbracht et all 120061 : 
iBeirao etal 1 l2008l:lGalliano et al.ll2008l) . 

Electron and ion interactions with dust and the implications 
of those interactions for the dust evolution and emission have 
alread y been d i scusse d in the litera t ure (e . g . iDraine & Salpeter 
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Jones et alJl996l) . In this work we extend this earlier work to the 
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case for PAHs, using our study of PAH evolutio n due to ion and 
electr on interactions in shock waves in the ISM dMicelotta et al.l 
I2009L hereafter MJT). Here we consider the fate of PAHs in 
the hot gas behind fast non-radiative shocks and in a hot gas in 
general. 
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The aim of this paper is to study the PAH stability against 
electron and ion collisions (H, He and C) in a thermal gas with 
temperature T in the range 10 3 - 10 8 K. 

The paper is organized as follows: Sect. 2 and Sect. 3 de- 
scribe the treatment of ion and electron interactions with PAHs, 
Sect. 4 illustrates the application to PAH processing in a hot gas 
and Sect. 5 presents our results on PAH destruction and lifetime. 
The astrophysical implications are discussed in Sect. 6 and our 
conclusions summarized in Sect. 7. 



2. Ion interaction with PAHs 

2.1. Electronic interaction 

The ion - PAH collision can be described in terms of 
two simultaneous pro cesses which can be treated separately 
(iLindhard et al.l fl963l) : nuclear stopping or elastic energy loss 
and electronic stopping or inelastic energy loss. 

The nuclear stopping consists of a binary collision between 
the incoming ion (projectile) and a single atom in the target ma- 
terial. A certain amount of energy will be transferred directly to 
the target atom, which will be ejected if the energy transferred 
is sufficient to overcome the threshold for atom removal. The 
physics of the nuclear interaction for a PAH target was presented 
in the companion paper (MJT). A summary of the theory is pro- 
vided in Sect. 2.2 and we here present the results of our calcula- 
tion. 

In this paper we focuse on the electronic stopping, which 
consists of the interaction between the projectile and the elec- 
trons of the target PAH, with a subsequent energy transfer to the 
whole molecule. The resulting electronic excitation energy will 
be transferred to the molecular vibrations of the PAH through ra- 
diationless processes (eg., interconversion & intramolecular vi- 
brational redistribution). The vibrationally excited molecule will 
decay through either (IR) photon emission or through fragmen- 
tation (i.e., H-atom or C2H, loss, where n = 0, 1, 2). 

No specific theory describes the energy transfer to a PAH via 
electro nic excitation, so we adop t the same approach deve loped 
by e.g. ISchlatholter et all (1 19991) and lHadiar et al.l (1200 ll) who 
modelled electronic interactions in fullerene, Ceo- To calculate 
the energy transferred to a PAH, we treat the large number of 
delocalized valence electrons in the molecule as an electron gas, 
where the inelastic energy loss of tra veling ions is due to long 
range coupling to electron-hole pairs dFerrellll 19791) . 

In the energy range we consider for this study, the energy 
transferred scales linearly with the velocity v of the incident ion 
and can be described in term of the stopping power S , which is 
widely used in the treatment of ion-solid collisions. The stopping 
power represents the energy loss per unit length and is defined as 
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where dT is the energy loss over the pathlength ds dSigmundl 
[1981 . The total energy loss can then be obtained integrating 
Eq.ffl 



y(r s ) vds 



(2) 



The friction coefficient y dPuska & Nieminenll 19831) depends on 

the density parameter r s = (|7r«o) ^ which is a function of the 
valence electron density «o- The choice of the valence electron 
density is a delicate issue for PAHs. In the case of fullerene, this 




Fig. 1. The coordinate system adopted to calculate the energy 
transferred to a PAH via electronic excitation by ion collisions 
and by impacting electrons. The molecule is modeled as a disk 
with radius R and thickness d. The trajectory of the incoming 
particle is identified by the angle §, while the angle a corre- 
sponds to the diagonal of the disk. 



is assu med to be the spherically sy mmetric jellium shell calcu- 
lated bv lPuska & Nieminenl ( 1993b. which can b e approximated 
by the following expression ( Hadiar et al.ll2.00lT> : 



n = 0.15 exp[-(r-6.6) 2 /2.7] 



(3) 



where r (in atomic units, a.uQ) is the distance from the fullerene 
center. The electron density decays outside the shell and toward 
the center in such a way that the major contribution comes from 
the region with 4 < r < 9. 

The similarity in n electronic structure and bonding allows us 
to apply this jellium model also to PAHs. However, the spherical 
geometry is clearly not appropriate for PAHs, which we model 
instead as a thick disk analogously to the distribution from Eq. [3] 



n = 0.15 exp [-(xf / 2.7] 



(4) 



where x is the coordinate along the thickness of the disk. This 
electron density peaks at the center of the molecule and vanishes 
outside, leading to a thickness d ~ 4.31 A. The radius R of the 
disk is given by the usual expression for the radius of a PAH: 
QpAH = R = 0.9 yf(Nc) A, where Nr . is the number of carbon 
atoms in the molecule (lQmontlll986l) . For a 50 C-atoms PAH, 
R = 6.36 A. 

To calculate the energy transferred from Eq. [2] we adopt the 
coordinate system shown in Fig. 1, where the pathlength through 
the PAH s is expressed as a function of the coordinate x and of 
the angle § between the axes of the molecule and the direction 
of the incoming ion. In this way for each trajectory given by § 
the corresponding energy transferred can be computed. We have 
d* -ds cos § and the electron density is then given by 



«oO, ft) = 0.15 exp [-(s cos ff) 2 / 2.7 ] 



(5) 



The density parameter can be rewritten as r s (s, 1?) = 
(^nna(s, 1?)) ' The friction coefficient y has been calculated 



1 Length (a.u.) = a = 0.529 A, Velocity (a.u.) = 0.2 V£(keV/amu), 
Energy (a.u.) = 27.21 16 eV 
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bv lPuska & Nieminenl (1 19831) for various projectile ions and r s 
values. It can be interpolated by the exponential function 



yOs) = To exp 



-(r s (s,0)-1.5) 



R 2 



(6) 



where To is the value of y when r s = 1.5. For hydrogen, he- 
lium and carbon the fit parameter R 2 equals 2.28, 0.88 and 0.90, 
respectively. 

The energy transferred is then given by the following equa- 
tion 



T e (ff) =27.2116 



XR/sin i? 
vy(r s )ds 
Rj sin i? 



(7) 



with y(r s ) from Eq.[6] The constant 27.21 16 converts the dimen- 
sions of T e from atomic units to eV. For a given velocity of the 
incident particle, the value of T e will be maximum for i? = n/2 
(s = 2R, the longest pathlength) and minimum for § — 0, n 
(s = d, the shortest pathlength). The deposited energy also in- 
creases with the pathlength across the molecule, and then will 
be higher for larger PAHs impacted under a glancing collision. 

The energy transferred determines the PAH dissociation 
probability upon electronic excitation, which is required to quan- 
tify the destruction induced by inelastic energy loss in the hot gas 
(see Sect. 4). 

2.2. Nuclear interaction above threshold 

In the present study we have to consider not only the electronic 
interaction described above, but the nuclear part of the ionic col- 
lision as well. The theory of nuclear interaction above threshold 
has been described in detail in MJT We summarize here for clar- 
ity the essential concepts and the equations which will be used 
in the following. 

For the nuclear interaction we consider only collisions above 
threshold, i.e. collisions able to transfer more than the minimum 
energy Tq required to remove a C-atom from the PAH. To calcu- 
late the PAH destruction due to nuclear interaction, we use the 
rate of collisions above threshold between PAHs and ions in a 
thermal gas, as given by Eq. 30 in MJT: 



J 

•JVi) 



R Bj = N c 0.5xinn F c vcr(v)f(v,T)dv 



(8) 



where /(v, T) is the Maxwellian velocity distribution function 
for ion, i, in the gas, Fq is the Coulombian correction factor, <x(v) 
is the nuclear interaction cross section per atom above threshold 
(Eq. 22 in MJT where E is the kinetic energy corresponding to 
v), and the factor 1/2 takes the averaging over the orientation of 
the disk into account. 

For To we adopt a value of 7.5 eV (see discussion in MJT, 
Sect. 2.2.1). The kinetic energy required for the incoming ion 
to transfer Tq is the critical energy Eo. The lower integration 
limit v'o in Eq. [8] is the critical velocity corresponding to Eo a , 
which is the minimum kinetic energy for the projectile to have 
the nuclear interaction cross section different from zero (cf. Sect. 
2.2 in MJT). 

3. Electron collisions with PAHs 

Fast electrons are abundant in a hot gas. Because of their low 
mass, they can reach very high velocities with respect to the 
ions, and hence high rates of potentially destructive collisions. 
We consider gas temperatures up to 10 8 K, corresponding to a 
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Analytical fit to the electron stopping power in solid 
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thermal electron energy of ~10 keV, well below the relativis- 
tic limit of ~500 keV. Under these low-energy conditions, with 
respect to the relativistic regime, elastic collisions between elec- 
trons and target nuclei are not effective. The energy transfer oc- 
curs through inelastic interactions with target electrons (as for 
electronic excitation by impacting ions), which lead to a collec- 
tive excitation of the molecule, followed eventually by dissocia- 
tion or relaxation through IR emission. 

The calculation of the energy transferred by such 'slow' elec- 
trons is in fact a delicate matter. The t heory developed under the 
first Born approximation dBefhd 19301) can be applied only to the 
most energetic electrons (around few keV) but is unsuitable for 
the rest of our range. At low energies (<10 keV), where the first 
Born approximation is no longer valid, the Mott elastic cross 
section must b e used instead of the conventional Rutherfor d 
cross section dMott & Massevl 119491: ICzvzewski et al . 1990). 
The Monte Carlo program CASINO dHovington et al.l 119971) 
computes the Mott cross sections in the simulation of electron 
interactions with various materials. Unfortunately the stopping 
power dE/dx is not included in the program outp ut. An em- 
pirical expression for dE/dx has been proposed by iJov & Luol 
(fl989l) . which nevertheless is reliable only down to 50 eV, while 
we are interested in the region between 10 and 50 eV as well. 

We decided thus to derive the electron stopping power from 
experimental results. Measurements of the electron energy loss 
in PAHs are not available in our energy range of interest, so we 
use the measurement of dE/dx in solid carb on for electrons with 
energy between 10 eV and 2 keV dJovll995i) . The data points are 
well fitted (to within few %) by the following function: 



S(E) = 



h log(l +aE) 
fE'i + bE d + cE e 



(9) 



where E is the electron energy (in keV). The values of the fitting 
parameters are reported in TableQ] S (E) has the same functional 
form as the ZBL reduced stopping cross section for nuclear in- 
teraction (cf. Sect. 2.1 in MJT). The datapoints and the fitting 
function are shown in the top panel of Fig. [2] 

The stopping power increases sharply at low electron ener- 
gies, reaches its maximum at ~0. 1 keV, and decreases smoothly 
afterwards. Between 0.01 and 0.1 keV a small variation of the 
energy of the incident electron will translate into a large change 
in the transferred energy per unit length. The shape of S (E) im- 
plies that only those electrons with energies that fall in a well 
defined window will efficiently transfer energy, while electrons 
below ~0.02 keV and above ~2 keV are not expected to con- 
tribute significantly to PAH excitation. 

Once the stopping power dE/dx is known, we can calculate 
the energy transferred by an electron of given energy when trav- 
elling through the PAH. We adopt the same configuration used 
for electronic interaction, shown in Fig.Q] The trajectory of the 



4 



E. R. Micelotta et al.: PAH processing in a hot gas 



03 

o 



Q. 
Q_ 
O 

CO 





Joy 1 995 ■ 






J 


\ Carbon: S(E) - 

\ 









i 



0.01 



0.1 



Electron energy (keV) 



10 



•< 
ClT 




200 



Electron energy (eV) 



Fig. 2. Top panel - Experimental measurement of dE/dx in solid 
carbon for electrons with energy between 10 eV and 2 keV, from 
iJovl ( 1 19951) . overlaid with the fitting function S(E) (solid line). 
Bottom panel - Integral F calculated numerically as a function 
of the energy of the incident electron, E. 

incoming electron is defined by its impact angle § and by the 
geometry of the molecule. 

Because the thickness of the PAH is non-negligible with re- 
spect to its radius, the stopping power is not constant along the 
electron path. To calculate the energy loss we thus follow the 
procedure described below. We have dE/dx = -S(E), then 



/"*--/ 



dE 

S(E) 



F(E) 



(10) 



Thus, x\ — xn = - [F{E\) _ F(Eq)] where x\ — x = is 
the maximum pathlength through the PAH. From inspection of 
Fig-El one can see that, if | tan(i?)| < tan(ff), = d/|cos#|, 
otherwise l(§) = 2R/\ sini?|. Eq is the initial energy of the im- 
pacting electron, E\ is the electron energy after having travelled 
the distance l{ff), which needs to be calculated so we can then 
determine the energy transferred to the PAH 



r e lec(#) -E -Ei 



(ID 



The integral F, calculated numerically as a function of E, is 
shown in the bottom panel of Fig. [2] We recognize that for low 
energies, F rises sharply, reflecting the small energy stopping 
power in this energy range (cf., top panel in Fig. [2J. For higher 
energies, F(E) rises slowly (and linearly) with increasing energy 
over the relevant energy range. We note that the initial rise de- 
pends strongly on the (uncertain) details of the stopping power 
at low energies. However, it is of no consequence in our deter- 
mination of the amount of energy deposited since we are only 



concerned with those collisions for which the energy deposition 
is in excess of the threshold energy (^ 10 eV) over a pathlength 
of =i 10 - 20 A. Or to phrase it differently, at low energies a 
much larger pathlength has to be traversed (than is relevant for 
PAHs) in order to transfer sufficient energy to cause fragmen- 
tation. Then, since we know Eq, we can calculate F(Eq), and 
F(Ei) is then given by F(Eq) - We also calculated £ as a 
function of F(E). We can then determine E\ from F{E\) and, 
finally, r e i ec from Eq. [TT] 

4. PAH destruction 

4.1. Dissociation probability 

Ion or electronic collisions (or UV photon absorption) can 
leave the molecule internally (electronically) excited with an en- 
ergy Te- Internal conversion transfers this energy to vibrational 
modes, and the molecule can then relax through dissociation or 
IR emission. These two processes are in competition with each 
other. To quantify the PAH destruction due to ion and electron 
collisions we need to determine the probability of dissociation 
rather than IR emission. 

In the microcanonical description of a PAH, the internal en- 
ergy, Te, is (approximately) related to the effective temperature 
of the system, r e ff, by the following equation 



^eff 



2000 



r e (eV) 
N c 



1 - 0.2 



£o(eV) 
r e (eV) 



(12) 



where Eq is the binding energy of the fragment (lTielensll2005K 
This effective temperature includes a correction for the finite 
size of the heat bath in the PAH. The temperature r e ff is de- 
fined through the unimolecular dissociation rate, k& ss > written in 
Arrhenius form 



^diss = ko (T e s) exp [- E Q /k T eS ] 



(13) 



where k - 8 .617x10 5 eV/K is the Boltzmann's constant (cf. 
lTielensll2005l) . 

Consider the competition between photon emission at a rate 
AjR (photons s _1 ) and dissociation at a rate of fcdiss (fragments 
s _1 ). For simplicity, we will assume that all photons have the 
same energy, As. The probability that the PAH will fragment 
between the nth and (n + l)th photon emission is given by 



<Pn = Pn+i ]~~[(1 - Pi) 

i=0 



(14) 



The (un-normalized) probability per step p-, is given by = 
kdiss{Ei)/kiR(Ej) and £; = (Te-JX As), with Te the initial internal 
energy, coincident with the energy transferred into the PAH. The 
total dissociation probability is then given by 



f(«max) = 2^ f" 



(15) 



n=0 



If we ignore the dependence of km on the temperatures 
(Eq.fT2l drop approximately by 7y7Vi = (1 - 0AAs/Te). The 
probability ratio is approximately: 



Pi 
Pi-i 



exp 



Eq 
kT 



0.4- 



(16) 



These equations become very difficult to solve in closed form. 
However, let us just assume that /?,- does not vary and is given 
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by pav Then we have the total un-normalized dissociation prob- 
ability 



20 



■P("max) = ("max + l)Pav = 



k exp [-E /k ± avj 

^IR/(«max + 1) 



(17) 



where we adopt &ir = 100 photons s 1 dJochims et al.ll 19941) . and 
the average temperature is chosen as the geometric mean 



(18) 



where Tq and T„ mia are the effective temperatures corresponding 
to Te (initial internal energy equal to the energy transferred) and 
(Te - « m ax x As) (internal energy after the emission of « max pho- 
tons). For the energy of the emitted IR photon we adopt the value 
Ae = 0.16 eV, corresponding to a typical CC mode. 

If the p/s were truly constant, then n max would be « max = 
{Te - Eq)/Ae. However, they do decrease. So, rather we take it 
to be when the probability per step has dropped by a factor 10. 
A direct comparison between the full evaluation and this simple 
approximation yields « max = 10, 20 and 40 for Nc = 50, 100 and 
200 respectively. The quantity « max scales with Nc because for 
a constant temperature (eg., required to get the dissociation to 
occur), the internal energy has to scale with Nc- As a result, the 
number of photons to be emitted also has to scale with Nc- 

The choice of the values to adopt for ko and Eq is a deli- 
cate matter. In the laboratory the dissociation of highly vibra- 
tionally excited PAHs is typically measured on timescales of 1- 
100 ps because either the molecules are collisionally de-excited 
by ambient gas or the molecules have left the measurement zone 
of the apparatus. In contrast, in the ISM, the competing relax- 
ation channel is through IR emission and occurs typically on a 
timescale of 1 s. As is always the case for reactions characterized 
by an Arrhenius law, a longer timescale implies that the internal 
excitation energy can be lower. This kinetic shift is well estab- 
lished experimentally and can amount to many eV. Moreover, 
only small PAHs (up to 24 C-atoms) have been measured in 
the laboratory and the derived rates have to be extrapolated to 
much larger (~ 50 C-atoms) PAHs that are astrophysically rele- 
vant. In an astrophysical context, the unimolecular dissociation 
of highly vibrationally excited PAHs - pumped by FUV pho tons 
- has been studied experimentally by Jochims et alj (fl994l) and 
further analyz ed bvlLe Page et al.l (120011) . Here, we will modify 
the analysis o flTielensI d2005l) for H-loss by UV pumped PAHs to 
determine the parameters for carbon loss. The dissociation rate - 
given by Eq.[T_3]- is governed by two factors, the pre-exponential 
factor ko and the energy Eq. The pre-exponential factor is given 
by 



kT eS 



exp 



1 + 



AS 
~R~ 



(19) 
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where A S is the entropy chan ge which we set equal to 10 cal 
mole 1 riLing&Lifshitzlll998h . We will ignore the weak tem- 
perature dependence of ko in the following and adopt 1.4 x 10 16 
s _1 . The param eter Eq can now be d etermined from a fit to the 
experiments bv lJochims et alj (1 19941) on small PAHs. These ex- 
periments measured the appearance energy £ ap - the internal en- s~ l in a Habing field). Hence, if the lifetime of the smallest PAHs 



Fig. 3. The appearance energy as a function of the number of C- 
atoms in the PAH. The red line provides a fit to the experimental 
data using Eq. [13] for an assumed pre-exponential factor ko = 
1.4 x 10 16 and a (fitted) Arrhenius energ y, E - 3.65 eV Th e 
data points are the experimental results of iJochims et al.l (Il994l) . 
The blue line is the appearance energy for ISM conditions (eg., 
at a rate of 1 s~'). 



apparent. Indeed, this kinetic shift can amount to tens of eV for 
relevant PAHs. 

The derived Arrhenius energy of 3.65 eV is small com- 
pared to the binding e nergy of a C2H2 group in a PAH (4.2 eV, 
iLing & Lifshitzlll998l) . This is a well known problem in statis- 
tical unimolecular dissociation theories (cf. iTielensI l2008h . We 
emphasize that these results show that a typical interstellar PAH 
with a size of 50 C-atoms would have a dissociation probabil- 
ity of ~ 1/2 after absorption of an FUV photon of ~12 eV (cf. 
Fig. |3j. Hence, PAHs would be rapidly lost in the ISM through 
photolysis. It seems that the experiments on small PAHs can- 
not be readily extrapolated to larger, astrophysically relevant 
PAHs. Possibly, this is because experimentally C2H2 loss has 
only been observed for very small catacondensed PAHs with a 
very open carbon skeleton (e.g., naphthalene, anthracene, and 
phenanthrene) which are likely much more prone to dissocia- 
tion than the astrophysically more relevant pericondensed PAHs. 
Indeed, the small pericondensed PAHs, pyrene and coronene 
did not show any d i ssociation on the exp erimental timescales 
(IJochims et alJI 1 994t One & Lifshitzll 1 998t> . 

Turning the problem around, we can determine the Arrhenius 
energy, Eq, as a function of the dissociation probability by adopt- 
ing an IR relaxation rate of 1 s -1 and an internal excitation en- 
ergy equal to a typical FUV photon energy (12 eV). The results 
for a 50 C-atom PAH are shown in Fig [4] If we adopt a lifetime, 
tpah, of 100 million years, a PAH in the diffuse ISM will have 
typically survived some ct uy N uy tpah = 2 x 10 6 /Vc UV photon 
absorptions (with <x uv = 7 x 10~ 18 cm 2 , N uv = 10 8 photons crrT 2 



ergy at which noticable dissociation of the PAH first occurred. 
The rate at which this happened was assumed to be 10 4 s _1 . 
Adopting this value, the appearance energy can be estimated 
from Eq.Qj] The results of our analysis are shown in Fig. [3] The 
internal energy required to dissociate a PAH depends strongly on 
the PAH size. Likewise the kinetic shift associated with the rel- 
evant timescale at which the experiment was performed is quite 



in the ISM (eg., with Nc - 50 C-atom) is set by photodissocia- 
tion of the C-skeleton, the probability for dissociation has to be 
5 x 10~ 7 corresponding to an Arrhenius energy of 4.6 eV (Fig. 
gj. We note that in a PDR the photon flux is higher (G ~ 10 5 ) 
while the lifetime (of the PDR) is smaller (t pdr ^ 3 x 10 4 yr), 
resulting in 6 x 10 7 UV photons absorbed over the PDR lifetime. 
Survival of PAHs in a PDR environment would therefore require 
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Fig. 4. The probability for dissociation of a 50 C-atom PAH ex- 
cited by 10 eV as a function of the Arrhenius energy, Eo- 

a somewhat larger £0 (or alternatively, only slightly larger PAHs 
could survive in such an environment). 

We note that the binding energy of a C2H2 group to small 
PAHs is estimated to be 4.2 eV and is probably somewhat larger 
for a 50 C-atom condensed PAH. Loss of pure carbon, on the 
other hand, requires an energy of 7.5 eV, close to the bind- 
ing energy of C to graphite. Loss of C2 from ful lerenes has a 
measured E of 9.5 + 0.1 eV dTomita et al.ll200lh . These latter 
two unimolecular dissociation channels are for all practical pur- 
poses closed under interstellar conditions. Finally, likely, H-loss 
will be the dominant destruction loss channel f or large PAHs 
(Eq = 3.3 eV l lead ing to rapid loss of all H's dLe Page et al.l 
l200lt lTielens|[2005h . The resulting pure C-skeleton may then 
isomerize to much more stable carbon clusters, in particular 
fullerenes, and t his may be the dominant 'loss' channel for inter- 
stellar PAHs (cf. lTielensl20"08l and references therein). It is clear 
that there are many uncertainties in the chemical destruction 
routes of interstellar PAHs and that these can only be addressed 
by dedicated experimental studies. For now, in our analysis of 
the unimolecular dissociation of PAHs - excited by electron or 
ion collisions - we will adopt Eq = 4.6 eV as a standard value. 
We will however also examine the effects of adopting £0 = 3.65 
eV (eg., p = 1/2) and E = 5.6 eV (eg., p = 3 x 10~ 12 ). 

4.2. Collision rate 

Once the dissociation probability is determined, we can calcu- 
late the destruction rate through electronic excitation following 
electron or ion collision. Adopting the configuration shown in 
Fig-H the destruction rate is given by 

fle = v X n H F c ^J dO (7g(0) P(v, 0) (20) 

with O = sin §d-&d(p, <p running from to 2n and # from to 
n/2. Eq.|20]can be re-written as 

Re = vxnnFc o- g (#)P(v,#) sin#d# (21) 
J#=o 

The term v is the velocity of the incident particle, x is the 
ion/electron relative abundance in the gas, and Fq is the 
Coulombian correction factor, taking into account the fact that 



both target and projectiles are charged, and thus the collision 
cross section could be enhanced or diminished depending on the 
charge sign. The electron coulombian factor is always equal to 
1 (to within 1%) because electrons have low mass. For a fixed 
temperature they reach higher velocities, with respect to the ions, 
and are thus less sensitive to the effects of the Coulombian field. 
The quantity # is the angle between the vertical axes of the PAH 
and the direction of the projectile. The term cr g is the geometrical 
cross section seen by an incident particle with direction defined 
by )?. The PAH is modelled as a thick disk with radius R and 
thickness d. The cross section is then given by 

o-g = nR 2 costf + 2Rd sin)? (22) 

which reduces to cr g = nR 2 for & — (face-on impact) and to 
cr g = 2Rd for 1? = n/2 (edge-on impact). The term P(v, §) rep- 
resents the total probability for dissociation upon electron colli- 
sions and electronic excitation, for a particle with velocity vpah 
and incoming direction §. It was calculated from Eq. [17] using 
the appropriated value of the energy transferred: r e i ec for elec- 
trons (Eq. ITTb and T e for electronic excitation (Eq. [7J 

We are considering a hot gas and therefore we are interested 
in the thermal collision rate, given by 

fl e .T= R e (v)f(v,T)dv (23) 

where vq e is the electron/ion velocity corresponding to the elec- 
tronic dissociation energy Eq. The temperature T is the same for 
both ions and electrons, but these latter will reach much larger 
velocities. From Eq. [2J] we then expect to find a significantly 
higher collision rate with respect to the ion case. 

5. Results 

5. 1 . PAH destruction in a hot gas 

To describe the destructive effects on PAHs of collisions with 
ions i (i = H, He, C) and electrons in a thermal gas we have to 
evaluate the rate constant for carbon atom loss. 

For all processes the rate constant / is defined by the ratio 
J = Rilrin, where Rj is the thermal collision rate appropriate 
for nuclear and electronic excitation and electron collisions (cf. 
Eq. [8] and l23l , and «h is the density of hydrogen nuclei. For 
electronic and electron interactions the rate must be multiplyied 
by a factor of 2, to take into account the fact that each collision 
leads to the loss of two carbon atoms. 

The electron, nuclear and electronic rate constants for three 
PAH sizes (Nc = 50, 100 and 200 C-atom) are shown in Fig. 
[5] as a function of the gas temperature. We assume for the nu- 
clear threshold energy To the value of 7.5 eV (MJT), and for the 
electronic dissociation energy the value 4.58 eV (Sect. 4.1). 

For the nuclear interaction, the threshold energy To is inde- 
pendent from the PAH size, so the three curves start at the same 
temperature (not shown in the plot). The small separation be- 
tween the curves is due to the fact that the PAH 'surface area' - 
and hence the rate constant - scales linearly with Nc, therefore 
J is higher for bigger PAHs. 

For nuclear (and electronic) interactions, the rate constants 
decrease from hydrogen to carbon because of the lower abun- 
dance of the heavier projectiles with respect to hydrogen (H : He 
: C = 1 : 0.1 : 10~ 4 ). On the other hand, the nuclear curves shift 
toward lower temperatures from lighter to heavier projectiles. 
This is the reflection of the fact that the critical energy of the 
particle, required to transfer the threshold energy To, decreases 
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Fig. 5. Nuclear (dashed lines), electronic (solid lines) and electron (dashed-dotted lines) rate constant for PAH carbon atom loss due 
to collisions with H, He, C and electrons in a thermal gas. The rate constants are calculated as a function of the gas temperature 
for three PAH sizes Afc = 50, 100, 200 C-atoms, assuming the nuclear threshold energy Tq = 7.5 eV and the electronic dissociation 
energy Eq = 4.58 eV. 
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Fig. 6. Analytical fits to the calculated rate constants shown in Fig. [5] The Tonic' curve is the fit to the sum of the nuclear and 
electronic rate constants, thus represents the total contribution from ion collisions to PAH destruction. 
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with increasing mass of the projectile itself. Then, a carbon atom 
with a temperature of ~ 4x 1 3 K is hot enough to transfer the en- 
ergy required for atom removal via nuclear interaction, whereas 
for hydrogen a temperature of at least 10 4 K is necessary. The 
almost-constant behaviour after the initial rise reflects the large 
maximum observed in the nuclear cross section (cf. Fig. 2 in 
MJT). 

The dissociation probability P(n max ) (Eq.fTTTi depends on the 
binding energy of the fragment Eq, on the PAH size Nc and on 
the energy transferred (through T- dv ), which in turns depends on 
the initial energy (velocity) of the projectile. For a thermal dis- 
tribution, this latter will be determined by the gas temperature T. 

The electronic rate constant curves for the different PAH 
sizes are well separated at the lowest temperatures. This reflects 
the fact that, for a fixed value of the transferred energy and of 
the electronic dissociation energy Eq, the dissociation probabil- 
ity decreases for increasing Nc because either more energy is 
required in the bond that has to be broken or because the energy 
is spread over more vibrational modes and hence the internal 
excitation temperature is lower. On the other hand, the more en- 
ergy that is deposited into the PAH, the higher is the dissociation 
probability. The energy transferred via electronic excitation (and 
then r av ) increases with the energy of the projectile up to a max- 
imum value, corresponding to an incident energy of 100 keV for 
H (and higher for more massive particles), and decreases beyond 
that for higher energies. The energy content of a thermal gas at 
T — 10 8 K is ~9 keV, thus, in the temperature range considered 
in this study, the energy transferred increases with temperature 
(energy) and hence the dissociation probability increases as well. 
This is the basis for the monotonic rise of the electronic rate con- 
stant. After the initial separation, the three curves converge, be- 
cause the rise in the transferred energy compensates the effect of 
increasing A^c- 

As discussed in Sect. 3, the energy transferred by impacting 
electrons rises sharply for energies in excess of 10 eV, peaks at ~ 
100 eV and decreases more slowly down to 10 keV. This results 
in a dissociation probability shaped as a step function: for 10 eV 
S) v S= 4 keV, P(« max ) jumps from values close to zero up to 1. 
These limiting energies apply to a 50 C-atom PAH; for A^c = 200 
the width of the step is smaller (100 eV ^ v ^ 1 keV), due to the 
fact that for a bigger PAH more energy has to be transferred for 
dissociation. This behaviour is reflected in the shape of the elec- 
tron rate constant, where a steep rise is followed by a maximum, 
which is emphasised by the logarithmic scale used for the plot. 
As expected, the electron rate constant overcomes the electronic 
one, except for the highest gas temperatures. This results from 
the fact that, for a given temperature, electrons can reach higher 
velocities with respect to the ions (cf. Eq.l2"Tlandl2"3"ll. 

To summarize, from Fig.[5]we can infer that, according to our 
model, the destruction process is dominated by nuclear interac- 
tion with helium at low temperatures (below ~3xl0 4 K), and by 
electron collisions above this value. Small PAHs are easier de- 
stroyed than big ones for temperatures below ~10 6 K, while the 
difference in the destruction level reduces significantly for hotter 
gas. 

The calculated rate constants, shown in Fig. [5] are well fitted 
by the function g(T) = 10 /(log(r)) , where f(T) is a polynomial of 
order 5 



f(T) = a + bT + cT 2 + dT 3 +eT 4 +fT 5 



(24) 



For each PAH size we provide the fit to the electron rate constant; 
for the ions, the fit is over the sum of the nuclear and electronic 
rate constants, in order to provide an estimate of the global con- 
tribution from ionic collisions. The fits are shown in Fig. [6] and 



the fitting parameters are reported in Table [2] To provide an ex- 
ample of the accuracy of our fitting procedure, Fig.|7]shows the 
comparison between the calculated rate constant and the corre- 
sponding analytical fit, for electrons and helium impacting on a 
50 C-atom PAH. The He fit is over the sum of the nuclear and 
electronic rate constants. The average fitting discrepancy is ~15 
% and the fits are therefore well within any uncertainties in the 
method. 

Fig. [8] shows the comparison between the carbon loss rate 
constants for a very big PAH, Nc = 1000, and a 50 C-atom 
molecule. We assume To = 7.5 eV and Eq = 4.6 eV. Because of 
the decrease of the dissociation probability when the PAH size 
increases, as expected, the electron and electronic rate constants 
are strongly suppressed, and both curves shift toward higher tem- 
peratures. Indeed, for such a big PAH, a much higher internal en- 
ergy is required to reach the internal temperature where dissocia- 
tion sets in. On the other hand, the nuclear rate constant increases 
linearly with the PAH size. As a result, for a 1000 carbon atoms 
PAH, the nuclear interaction is the dominant (and efficient) de- 
struction mechanism up to T ~ 2xl0 7 K. In conclusion, elec- 
trons are responsible for the destruction of small/medium size 
PAH, while for big molecules this role is taken by ions. 

As mentioned at the end of Sect. 4. 1 , we examined the effects 
of adopting different values for the Arrhenius energy, Eq — 3.65 
and 5.6 eV, lower and higher, respectively, than our canonical 
value 4.6 eV. The results are shown in Fig. [9] The dissociation 
probability decreases for increasing Eq, because more energy is 
required in the bond that has to be broken. Hence, as expected, 
both the electron and electronic rate constants decrease in ab- 
solute value and shift toward highest temperatures. In particular 
the electronic thermal shift is very pronounced, indicating how 
sensitive this process is with respect to the assumed Eq. A varia- 
tion in the adopted electronic excitation energy translates into a 
significant variation of the rate constant, reemphasizing the im- 
portance of experimental studies of this critical energy. 

5.2. PAH lifetime 

Under the effect of electron and ion collisions in a hot gas, the 
number of carbon atoms in a PAH molecule varies with time. 
After a time t, this number is 

AW0 = AW0)exp[-f/T ] (25) 

and the number of carbon atoms ejected from this PAH is 

F L (0 = (l-exp[-//r ]) (26) 

The quantity to is the time constant appropriate for electron, nu- 
clear and electronic interaction, given by 



To 



N c 
(2)Rj 



Jn 



(27) 



H/e 



where Rj and J are the thermal collision rate and the rate con- 
stant for electrons, nuclear and electronic interactions respec- 
tively, and «H/ e is the hydrogen/electron density. For electrons 
and electronic excitation, the rate must be multipliyed by a factor 
of 2, because each interaction leads to the removal of two car- 
bon atoms from the PAH. The nuclear rate scales linearly with 
Nc, hence the corresponding fractional carbon atom loss F^ is 
independent of PAH size. 

For any given incoming ion and fixed PAH size Nc, F^ is 
univocally determined by the hydrogen/electron density rtn/e an d 
the gas temperature T. We assume that a PAH is destroyed after 
the ejection of 2/3 of the carbon atoms initially present in the 
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Fig. 7. Calculated rate constants for electrons and helium impacting on a 50 C-atom PAH, overlaid are the corresponding analytical 
fits. The He fit is for the sum of the nuclear and electronic rate constants. The average discrepancy is ~15 %. 

Table 2. Analytical fit parameters for the PAH carbon atom loss rate constant, calculated for electron and ion collisions. 





Fitting 


function 




f(T) = 


a + bT + c T 2 


+ d T 3 + e T 4 + f T 5 










a 


b 


c 


d 


e 


/ 






Electrons 


2136.83 


1632.17 


-499.822 


16.4341 


-5.82964 


0.177174 






H 


-1896.69 


1480.8 


-462.733 


71.8957 


-5.54719 


0.169996 


N c 


= 50 


He 


-971.448 


770.259 


-245.561 


38.8995 


-3.05787 


0.0954303 






C 


-704.392 


551.313 


-175.063 


27.6506 


-2.16871 


0.0675643 






Electrons 


-2255.38 


1681.45 


-503.451 


75.4339 


-5.64956 


0.168959 






H 


-1645.64 


1257.06 


-384.309 


58.4103 


-4.41007 


0.132356 


N c 


= 100 


He 


-945.901 


747.921 


-237.984 


37.6808 


-2.96613 


0.0928852 






C 


-711.244 


558.48 


-177.983 


28.2547 


-2.23145 


0.0701374 






Electrons 


-2234.37 


1597.71 


-459.647 


66.332 


-4.79841 


0.139019 






H 


-1473.64 


1109.01 


-334.292 


50.1417 


-3.74133 


0.111164 


N c 


= 200 


He 


-963.188 


765.639 


-245.054 


39.0745 


-3.10123 


0.0980047 






C 


-738.791 


584.928 


-187.884 


30.0819 


-2.39748 


0.0760639 



molecule. This occurs after the time tq (Eq.l27li which we adopt 
as the PAH lifetime against electron and ion bombardment in a 
gas with given density and temperature. 

Table [3] summarizes relevant data and PAH lifetime, To, for 
four objects characterized by warm-to-hot gas, X-ray emission, 
and (bright) IR emission features. The lifetime has been calcu- 
lated for PAHs with Nc = 200. To(ref) has been evaluated adopt- 
ing our reference values for the interaction parameters, Eq = 4.6 
eV and To = 7.5 eV, To(min) corresponds to the minimum values 
Eq = 3.65 eV and Tq = 4.6 eV, and ro(max) to the maximum 
values Eq = 5.6 eV and Tq = 15 eV. T bj ec t is the lifetime of the 
object. 

Clearly, PAHs or related larger species can survive in these 
environments. Fig. [TOl shows the fractional C-atom loss, due to 
electron and ion (H + He + C) collisions, for two widely differ- 
ent objects: the Orion Nebula (M 42) and the M 82 galaxy (cf. 
Table [3). The famous Orion Nebula is an Hn region with high 
density («h = 10 4 cm -3 ) and low temperature (T = 7000 K) gas, 
while M 82 is a starburst galaxy, which shows outside the galac- 
tic plane, a spectacular bipolar outflow of hot and tenuous gas 
(« H = 0.013 cm" 3 , T = 5.8xl0 6 K). 



In M 82, PAHs are completely destroyed by electrons, even 
for the larger PAH, before electronic and nuclear contributions 
start to be relevant. The electron and electronic fractional losses 
decrease with PAH size, thus bigger molecules can survive 
longer, while the nuclear loss is independent of iVc- The destruc- 
tion timescale is very short: after one thousand years the PAHs 
should have completely disappeared. 

In Orion the situation is reversed. At the temperature con- 
sidered for this object, electrons and electronic excitation do not 
contribute to PAH erosion (cf. Fig. [5). The damage is caused by 
nuclear interaction due to He collisions, with a marginal contri- 
bution from carbon because of low abundance, and the timescale 
is much larger: only after 10 million years the PAH destruction 
becomes relevant. Of course, we have not evaluated the destruc- 
tion of PAHs by H-ionizing photons in the Orion Hn region, 
which is expected to be very important. 

The young (^ 2500 yr) supernova remnant, N132D, in 
the Large Magellanic Clouds has been s tudied in detail at 
IR, optical, UV, and X-ray wavelengths dMorse et alJ 1 1995b 
iTappe etal]|2006l) . A Spitzer/IRS spectrum of the Southern rim 
shows evidence for the 15-20 fim plateau - often attributed 
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Fig. 8. Carbon atom loss rate constant for electrons and hydrogen impacting against a 1000 C-atom PAH. The rate constants for Nc 
= 50 are shown for comparison. 




Temperature (K) 

Fig. 9. Comparison between carbon atom loss rate constants calculated assuming three values for the electronic dissociation energy, 
Eq\ 3.65, 4.58 and 5.6 eV. The curves refer to electrons and hydrogen impacting against a 100 C-atom PAH. The nuclear rate 
constant is calculated assuming Tq = 7.5 eV. 



to large PAHs or PAH clusters jVan Kerckhoven etal l l2000t 
iPeeters et al.ll2004l) - and, tentat ively, weak PAH emission fea- 
tures n ear 6.2, 7.7, and 11.2 fjm dTappe et alJl2006l) . lTappe et alJ 
(120061) attribute these features to emission from large (~ 4000 
C-atom) PAHs either just swept up by the blast wave and not 
yet completely destroyed by the shock or in the radiative precur- 
sor of the shock. We calculate a lifetime of small (50 - 200 C- 
atom) PAHs in the relatively dense, hot gas of this young super- 
nova remnant of ~ 4 months (Table [3). In contrast, we estimate 
a lifetime of 150 yr for 4000 C-atom species (cf. Eq. [27] and 
Fig. [8}. It is clear that the PAH-grain size distribution will be 
strongly affected in this environment. Given advection of fresh 
material into the shocked hot gas, the observations are in reason- 
able agreement with our model expectations. We note that the 
observed shift to larger PAH sizes - so evident in the observa- 
tions - implies that the emitting species are not associated with 



the precursor but are actualy tracing the postshock gas and are 
therefore likely heated through collisions with the hot electrons. 

This discussion clearly shows that in a low density - high 
temperature gas, small PAHs are rapidly and completely de- 
stroyed. Survival of small PAHs in such environnement (and so 
the possibility to be detected) requires a protective environment 
and/or an efficient reformation mechanism. 

6. Discussion 

6.1. X-ray absorption 

As our calculations show, PAHs are rapidly destroyed in the hot 
gas associated with stellar winds and supernova explosions. Any 
PAHs observed near such regions have to be isolated from this 
hot gas and are presumably in cold gas entrained in these stellar 
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Table 3. Physical properties and PAH lifetime, To, in objects showing PAH emission characteristics. 



Object 




T 


Reference 


T (ref) 


T (min) 


To(max) 


^"object 


Destruction 




(cm- 3 ) 


K 




(yr) 


(yr) 


(yr) 


(yr) 


agent 


M82 


0.013 


5.8xl0 6 


(a) 


162 


43 


1302 


~2xl0 7 


electrons 


M17 


0.3 


7xl0 6 


(b) 


8.2 


2 


71 


~lxl0 6 


electrons 


N132D 


10 


8xl0 6 


(c) 


0.3 


0.07 


2.5 


~3xl0 3 


electrons 


Orion 


10 4 


7xl0 3 


(d) 


1.3xl0 7 


1.4xl0 5 


7.5x10" 


-lxlO 6 


ions 



(a): lRanalli etal] d2008h . (b): lTownslev et alJd2003l) . (c): lHwang et all jl993h . (d): lTielensl d2005h 



and galactic winds. However, such PAHs would still be exposed 
to en ergetic X-ray photons and these can be very destructive as 
well dVoitlll99llBoechat-Robertv et al.ll2008l) . 

The photon absorption cross section of PAHs shows strong 
peaks at about 6 and 17.5 eV associated with transitions involv- 
ing the n and cr electronic manifolds and then anothe r broad peak 
around 286 eV due to carbon K-shell transitio ns dKeller et al.l 
Il992t Ide Souza etal] 120021: iRegier et al.ll2007l) . Each of these 
peaks can show various subpeaks due to electronic and vi- 
brational structure. Here, we focus on the high energy peak. 
Because we are only concerned with the Planck-averaged cross 
section, all the fine detail will be washed out and we have elected 
to evaluate the X-ray absorption rate adopting the measured ab- 
solut e cross section of grap hitic carbon taken from the NIST data 
base (IChantlerll995ll2000l) . Fig.QT]shows the (photon) averaged 
cross sections as a function of the black body temperature. As 
this figure illustrates, even the strong absorption edge due to the 
carbon K-shell is washed out by this averaging process, justi- 
fying our neglect of the fine detail in the absorption cross sec- 
tions of individual PAHs. Our averaged cross section is also in 
good a greement with the recent study of iBoechat-Robertv et al.l 
(120081) for benzene. 

The photon absorption rate, Rx, of a PAH exposed to an X- 
ray photon flux, Nx is then given by, 



Rx = cr x (7x) JV X 



(28) 



with cr x (Tx) the Planck average absorption cross section at tem- 
perature, Tx- Consider now PAHs embedded in cold gas en- 
trained in the stellar wind of the ionizing stars of the Hn re- 
gion M 17. The observed X-ray luminos ity is 4 x 10 33 erg s' 1 
and the X-ray temperature is 7 x 10 6 K dTownslev et alJl2003h . 
With a projected emitting area of 42 pc 2 this then implies an av- 
erage X-ray photon intensity of ^ 2.6 x 10 3 photons cirT 2 s -1 . 
With an average cross section of cr x ^ 9 x 10~ 20 cm 2 /C-atom, 
we calculate a photon absorption rate of 7.2 x 10~ 3 photons/C- 
atom/Myr. After X-ray photon absorption, a PAH will typically 
lose a few C-atoms and, hence, the estimated lifetime of PAHs 
entrained in the stellar winds from M 17 is ~ 50 Myr. This is 
in excess of the stellar lifetime and the overall life time of the 
region. PAHs embedded in cold gas entrained in the galactic su- 
perwind of M 82 are exposed to somewhat harsher conditions. 
For M 82, at a distance of 3.6 Mpc, the observed soft (0.5-2.0 
keV) X-ray flux is 10~ n erg cuT 2 s and the temperature is 
7xl0 6 K. Evaluating the X-ray flux at a distance of 3 kpc from 
the nuclear starburst region, we calculate a photon absorption 
rate of 1.8 x 10~ 2 photons/C-atom/Myr. Hence, in this environ- 
ment, PAHs would be destroyed on a timescale of ~ 20 Myr, 
comparable to or larger than the starburst lifetime. It is clear that 
PAHs in the hot gas are predominantly destroyed by collisions 



with the hot gas and that X-ray photon absorption plays little 
role. 

Finally, we note that here we have assumed that every X- 
ray photon absorption will lead to photodissociation. Actually, 
H-shell electron ejection will be rapidly followed by the Auger 
effect filling the K-shell again and simultaneously ejecting a sec- 
ond electron. This process may leave the PAH internally excited 
with typically 15-35 eV0. This energy is well above the mea- 
sured appearance energy of PAH fragmentation for small PAHs 
under radiative cooling conditions (cf., Sect. 4.1) but larger 
PAHs may survive. Indeed, for an Arrhenius energy of 4.6 eV, a 
50 C-atom PAH requires about 24 eV of internal energy to dis- 
sociate with a probability of 1/2 (cf., Sect. 4. 1). Hence, not every 
X-ray photon absorption will lead to fragmentation, particularly 
for large PAHs. Further experiments on larger PAHs are required 
to settle this issue. 



6.2. PAHs as tracers of cold entrained gas 

The results in the previous subsection illustrate that PAHs can 
survive in the surface regions of cold gas clouds embedded in 
a hot gas for tens to hundreds of million years. The expansion 
of a stellar or galactic wind in an inhomogeneous environment 
will in a natural way lead to the entrainmen t of cold gas clumps 
dSmith et al.lll984tlChevalier & Clegd[l985l:lHartauist & Dvsoi ' 
1994t IStrickland & Stevensl 1200 
The material in these clumps will 



1 1993b 



ISuchkov et al.. 
120051) 



iMarcolini et al.l 

only be mildly shocked - and PAH destruction by shocks is 
not expected as long as the shock velocity is less than about 
100 km s _1 (MJT) - and can be transported with the wind while 
losing its identity only slowly, mainly through evaporation into 
the hot gas driven by thermal conduction. Multiwavelength 
observations of such winds support this global view and much 
of the soft X-ray and Ovi absorbing and emitting gas is now 
thought to originate from the interaction int erfaces between the 
tenuous hot gas and the co ld dense clumps ( IChevalier & Cleggl 
n^lHlckman et al.ll2002h . 

A typical 0.5 keV X-ray photon will penetrate a column of 
about 10 22 H-atoms/cm 2 and hence the 'PAH-bright' surfaces 
of such PDRs - which are typically only a column of 5 x 10 21 
H-atoms/cm 2 thick - will be slowly depleted of PAHs on the 
timescale calculated in section 16.11 However, at the same time, 
entrained clumps will slowly evaporate due to thermal conduc- 



2 The resulting PAH dication will quickly recombine with electrons 
and Coulomb explosion is of no consequence for such large PAHs. 
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Fig. 10. Nuclear and electronic fractional C-atom loss as a func- 
tion of time in the M 82 galaxy (top) and in the Orion Nebula 
(bottom), calculated for three PAH sizes N c = 50, 100, 200 C- 
atoms. The fractional carbon loss is defined as the number of 
carbon atoms ejected from the PAH after a given time divided by 
Nc, and represents a direct indicator of the level of destruction 
of the PAH. Each curve shows the total contribution from colli- 
sions with H, He and C. In the nuclear case the fractional loss 
is independent from Nc, so the three curves coincide, whereas 
for electrons and electronic interaction the fractional loss varies 
with the PAH size. In M 82 the PAH destruction is dominated 
by electron collisions and occurs on a very short timescale (~ 
thousand years). In Orion instead, because of the low gas tem- 
perature no electron and electronic contribution are possible (the 
rate constant is close to zero, see Fig. O, and hence only nuclear 
interaction is responsible for the PAH erosion, which however 
takes place on a long timescale (some 10 7 years). 



tion. Considering classical thermal conduction, the mass loss of 
a cloud due to evaporation into the hot gas is given by 



M 



(R \ I T \ 5 ^ 2 
z— I Moyr 1 (29) 
lOpcj \6x 10 6 K/ oy ' 



with T the temperature o f the hot gas and R the radius of the 
clump fcf. lTielens! 120051) . Consider the galactic wind of M 82 



o 

Si 




rm 



Fig. 11. Planck averaged cross section of PAHs as a function of 
temperature. The cross section is on a per C-atom basis. The 
strong carbon K shell band edge at 286 eV (2 x 10 6 K) has 
been largely washed out by the averaging process. 



and a typical structure with a size of 10 pc. Assuming thermal 
pressure equilibrium between the clump and the hot gas (eg., ig- 
noring ram pressure confinement which may actually be a factor 
of 10 larger), we estimate a density of 800 cirT 3 for a typical 
PDR temperature of 100 K for the galactic wind region of M 82. 
The lifetime of such a clump is then, 



m ( »pdr \( R \" / 6xl0 6 K \ 5/2 

Tdump * 40 (low ) 1 To^ ) [-^—j Myr ' (30) 

Hence, 'fresh' material will continuously be advected to the 
clump surface where ambient FUV photons can excite it. Thus, 
we conclude that such clumps can survive for a long time and 
that, proviso illumination by FUV photons, PAHs will form an 
excellent 'dye' for tracing the presence of cold entrained mate- 
rial; or more specifically, for the exposed surfaces of photodis- 
sociation regions. 



6.3. Comparison with previous studies 

iDwek et aD (1 19961) determines a sputtering rate for dust in a hot 
dust gas (T >10 6 K) that can be used to derive the lifetime, r of 
a grain of a given radius in a hot gas of a given density, i.e. 



10 3 / a 
«h \lA/ yi 



(31) 



where a is the grain size and «h is the gas density. Clearly in this 
situation it is the smallest grain that suffer the fastest destruc- 
tion. Hence, a treatment of the small grains, including PAHs, 
really needs to take into account the detailed physics of the in- 
teractions. For illustrative purposes, we can use Eq.[3T|to show 
that grains with sizes typical of PAHs (a as 5 A) should have life- 
times, in a tenuous hot gas («h = 1 cm -3 ), of the order of at least 
a few thousand years. We can see that this timescale is much 
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longer than the typical PAH lifetimes that we derive here, i.e. 
few years for a 50 carbon atom PAH in a gas with «h = 1 cm -3 . 
The timescale determination arrived at using Eq.[3T|assumes that 
a PAH just behaves as a small grain. However, in the case of a 
grain, the sputtering yield is usually much less that unity because 
atoms displaced by knock-on collisions in the solid can remain 
in the solid. The very much shorter lifetimes for the PAHs that 
we find can be ascribed to the fact that the transferred energy 
in incident electron and ion interactions always leads to the loss 
of carbon atoms from the PAH when that energy exceeds the 
relevant binding energy. In fact the equivalent PAH 'sputtering 
yield' can be greater than unity because multiple atom ejection 
is possible. Indeed we assume here that the ejected species in the 
case of electronic interactions are C2 units. 

The evolution of the dust and PAH thermal emission arising 
from a hot gas is affected by the chan ging dust size distribution. 
As pointed out bv lDwek et al.l d 1 996b . it is the short wavelength 
emission, coming from the smallest grains, that is most affected 
by dust destruction. Our work now indicates that the destruction 
timescale is much shorter than that predicted by Dwek et al. be- 
cause the erosion of the smallest grains and PAHs is enhanced 
by as much as three orders of magnitude com pared to the e arlier 
work. Thus, the gas cooling rates derived by iDwekl d 1987b will 
need to be reduced if a significant fraction of the dust mass is in 
the form of small grains and PAHs. 

6.4. C2 groups loss 

The primary channel for PAH erosion in a hot gas is the progres- 
sive loss of C2H, (« = 0, 1,2) units following incident electron 
and ion excitation of the molecule. We note that this can only 
really occur from the periphery of the PAH (see the lower part of 
Fig. 1 1 in MJT for an illustration of this type of erosion). In this 
case we can see that a coherent aromatic structure will tend to 
be preserved, which is probably not the case where the inertial 
sputtering of C atoms is dominant. Thus, we conclude that PAH 
erosion in a hot gas will tend to preserve the aromatic structure 
throughout the destruction process. 

As emphasized in MJT, PAH destruction may start with com- 
plete H-loss followed by isomerization to much more stable pure 
carbon clusters such as fullerenes (cf. iTiele ns 2008). While we 
have not assessed this point, we expect that such species will be 
more stable than PAHs in a hot gas. 

7. Conclusions 

We have extensively studied the stability of PAHs against elec- 
tron and ion collisions (H, He and C) in a hot gas, such as the 
gas behind fast non-radiative shocks and in galactic outflows. 
Collisions can lead to carbon atom loss, with a consequent dis- 
ruption and destruction of the molecule. We consider the case of 
a thermal gas with temperature T in the range 10 3 - 10 8 K. 

An ionic collision consists of two simultaneous processes 
which can be treated separately: a binary collision between the 
projectile ion and a single atom in the target (nuclear interaction) 
and energy loss to the electron cloud of the molecule (electronic 
interaction). 

For the nuclear interaction, as described in MJT, we have 
modified the existing theory in order to treat collisions able to 
transfer energy above a specific threshold To- We adopt To = 7.5 
eV as a reasonable value, but experimental determinations of this 
quantity are necessary. 

The electronic energy transfer has been described in term of 
the stopping power of an electron gas with appropriate electron 



density (jellium approximation). For electron collisions, we de- 
rived an analytical expression for the measured electron stopping 
power in graphite and applied this to the case of PAHs. 

The dissociation probability for a PAH excited by electronic 
interactions and electron collisions, has been derived using the 
theory of unimolecular reactions. The parameter Eq, which gov- 
erns the dissociation probability, is not well constrained. We 
adopt a value of 4.6 eV consistent with extrapolations to inter- 
stellar conditions but better determinations, relevant to the astro- 
physical situation, are needed. 

The PAH destruction process is dominated by electron col- 
lisions for gas temperatures above ~3xl0 4 K, and by nuclear 
interaction with helium below this value. Small PAHs are more 
easily destroyed than larger ones below ~10 6 K, while the dif- 
ference reduces significantly for a hotter gas. For a 1000 C-atom 
PAH, nuclear interactions are the dominant destruction process. 

In a hot and tenuous gas (e.g. M 82 galactic outflows), PAHs 
with sizes between 50 and 200 C-atom are destroyed by electron 
collisions in few thousand years. In denser and colder regions 
(e.g. Orion), PAHs can survive for some 10 7 yr before being de- 
stroyed by nuclear interaction processes. 

X-ray photon absorption plays little role in PAH destruction 
in the hot gas associated with stellar winds and supernova explo- 
sions, with respect to electron collisions. The PAH destruction 
process in a hot gas is therefore dominated by electron colli- 
sions. Any PAHs observed near such regions have to be isolated 
from this hot gas and are presumably in a cooler PDR-type gas 
entrained in the stellar and galactic winds. In this sense, PAHs 
represent a good tracer for the presence of entrained denser ma- 
terial. 

Our calculated PAH lifetime in a hot tenuous gas (T ~ 10 6 K, 
«h = 1 cm -3 ), is much shorter than the lifetime of an equivalent 
dust grain of roughly the same size (a « 5A). Thus, might then 
imply that the destructive effects of ion and electron collisions 
with very small grains have previously been underestimated. The 
enhanced erosion of the smallest grain and PAHs implies that the 
gas cooling rates for a hot gas (T > 10 6 K) - which depend on 
the grain/PAH size distribution - may need to be reduced if a 
significant fraction of the dust mass is locked in small grains and 
PAHs. 

PAH erosion in a hot gas occurs mainly through the ejection 
of C2 groups following electron collisions and electronic excita- 
tion. The C2 loss occurs at the periphery of the molecule, thus 
the aromatic structure will tend to be preserved throughout the 
destruction process. 
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